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Abstract 

The paper is devoted to the systematic derivation of nonperturba- 
tive thermodynamics of quark-gluon plasma, on the basis of the back- 
ground perturbation theory. Vacuum background fields enter only in 
the form of field correlators, which are known from lattice and ana- 
lytic calculations. In the lowest order in a s the purely nonperturbative 
sQGP thermodynamics is expressed through single quark and gluon 
lines (single line approximation) which are interacting nonperturba- 
tively with vacuum fields and with other lines. Nonperturbative EOS 
is obtained in terms (of the absolute value) of fundamental (adjoint) 
Polyakov loop Lf un M ad j\ and spatial string tension a s (T). In the lowest 
approximation the pressure for quarks (gluons) has a simple factorized 
form Pqtg) = PsBLf un d{adj)i where Li describe the action of vacuum 
colorelectric fields on particle trajectory. 

1 Introduction 

The study of the deconfined phase of QCD has a long history and was based 
mainly on the idea of weakly interacting gas of quarks and gluons, where 
perturbative and some nonperturbative methods can be used, (see a book in 
p] and reviews in [2]). 
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Recently, however, it was realized that there is a strong interaction in the 
QCD vacuum above the deconfmement temperature T c , as witnessed by the 
RHIC experiments, (see [3] for reviews) and lattice calculations (see e.g. [I] 
for a rewiew and further references). 

Independently of these developments, another approach was suggested in 
QCD almost two decades ago ([5], see [HI E| for a reviews), where the Nonper- 
turbative (NP) dynamics, instead of the perturbative one, was at the basis 
of QCD, and systematic formalism was developed both for T = ([HI [9]) 
and for T > ([10]- [12]). According to this approach, and in agreement 
with existing lattice data, the new picture of the deconfmed phase of QCD 
was suggested, where NP dynamics plays the dominant role. The latter is 
described by simplest (Gaussian) Field Correlators (FC) of color electric and 
color magnetic fields, D (x), D H (x), Df (x) and D^(x) which can be found 
from lattice measurements [131 EH] or from analytic calculations [15] . Here 
correlators D E (x), D H {x) are of pure NP origin, while D E ,H contain all per- 
turbative diagrams [16J. The phase transition was predicted to be due to the 
disappearance of D E (x) , which ensures colorelectric confinement, and the 
transition temperature T c was calculated (in terms of gluonic condensate) 
in good agreement with lattice data [17]. As a result above T c the NP dy- 
namics was predicted with colormagnetic confinement due to D H (x) and this 
prediction was supported later by detailed lattice measurements [T^ [T8]. 

This NP dynamics of (color)magnetic confinement plus the NP compo- 
nent of Df (which survives the phase transition) suggest a picture which is 
much different from the popular perturbative approach. In particular, the 
strong interaction in systems of quarks and gluons, leading to the formation 
of e.g. QQ bound state was predicted 15 years ago [19], and the role of the 
spatial string tension a s {T) for T >T C was stressed in [TUIIITJ, and supported 
by lattice data [EH El ED]. 

Recently, as a new challenge to the perturbative picture of quark-gluon 
plasma, the Debye masses have been computed through a s (T) [21], in good 
agreement with lattice data |22j . 

An important dynamical role in the formalism of field correlators is played 
by Polyakov loops, and the latter have been computed in terms of correlators 
and the property of Casimir scaling for them was theoretically established 
in [23]. Recently the Casimir scaling for five lowest representations SU(3) 
was accurately checked on the lattice [24J, thus supporting the quantitative 
accuracy of the method also for T > T c . At this point it is clear that a 
systematic NP formalism for all phenomena in T > QCD is needed, which 
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should be a basis for calculation of EOS, Polyakov loops, susceptibilities etc,. 
In the present paper, we undertake this task, giving a detailed derivation 
of the expressions for generating functions, pressure and internal energies 
of quarks and gluons, where perturbative contributions are considered as a 
series of Background Perturbation Theory (BPTh) and NP effects due to NP 
background can be taken fully in each term of 0(a"). 

The present paper is concerned with the lowest order 0{a° s ) of BPTh 
and gives the pressure of quarks and gluons in the form, where the leading 
contribution is due to single quark and gluon line interacting with the NP 
background (Single Line Approximation (SLA)), while line-line NP interac- 
tion is a correction. In SLA the pressure can be written in a simple factorized 
form Eqs. (66), (67), where NP factors are of the absolute value Polyakov 
loops, while colormagnetic background fields enter in the line-line correction. 

Our main region of interest is the low temperature region, T c < T < 2T C , 
where nonperturbative effects are most important, while at higher tempera- 
tures one expects possible connections of this formalism with the well devel- 
oped perturbation theory now known to 0(g 6 ) (see [25] and refs. therein). 

Two limitations are present in this paper: we do not consider chiral ef- 
fects, and possible accumulation of quark and gluon bound states near T c 

The plan of the paper is as follows. In section 2 the basics of BPTh and 
path integrals is given and free quark-gluon gas result is derived as an illus- 
tration. In section 3 NP interaction of quarks and gluons with background 
is expressed in the form of Polyakov loops (for colorelectric fields). Polyakov 
loops are studied in section 4 and EOS is summarized in section 5. Section 
6 is devoted to line-line interaction. Results and discussion are given in the 
concluding section. 

2 Background Perturbation Theory at T > 
2.1 Basic equations 

We start with standard formulae of the background field formalism [8j [9] 
generalized to the case of nonzero temperature in [10]- [12]. We assume that 
the gluonic field can be split into the background field and the (valence 
gluon) quantum field 

A fl = B ll + a M , (1) 
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both satisfying the periodic boundary conditions (PBC) 



B^{z 4 , Zi) = B^(z 4 + n(3,Zi); a M (z 4 , z { ) = a^(z 4 + n(3, z { ), (2) 
where n is an integer and f3 = 1/T. The partition function can be written as 

Z(V,T) = (Z(B)) B (3) 

with 

Z(B) =Nj Ztyexp (- J* dr J d 3 xL tot (x } r^j (4) 
and where <fi denotes all set of fields a M , L tot is defined in [10]- [12] and 



for the convenience of the reader in Appendix 1. N is a normalization con- 
stant. Furthermore, in Eq. (jSJ) ( )b means averaging over (nonperturbative) 
background fields B^. The precise form of this averaging is not needed for 
our purpose. 

Integration over the ghost and gluon degrees of freedom in Eq. yields 
the same answer as for the case T = [19J, but where now all fields are 
subject to the periodic boundary conditions ([^.Disregarding for simplicity 
quark terms and source terms in L tot in (j4[), one obtains 

Z{B) = N'(det W(5))-f [deti-D^D^B + a))U . 



Sint (a — jjj 



eg L , , '—-Jl 

I 



1 + E ""^",, \^v(-\jGj) , (5) 

1=1 L - I J^D^F^B) 

where the valence gluon Green's function G = W~ x is G^ = {D\ ■ 1 + 
2igF flv )~ 1 , and the tilde sign here and below refers to the operators in the 
adjoint representation, e.g. F h ^ v = iF^ v f abc . We can consider strong back- 
ground fields, so that gB^ is large (as compared to Aq CD ), while a s = g 2 /An 
in that strong background is small at all distances. Moreover, it was shown 
that a s is frozen at large distances[9]. In this case Eq. is a perturbative 
sum in powers of g n , arising from the expansion in {ga^) n . 

In what follows we shall discuss the Feynman graphs for the thermody- 
namic potential F(T,fi), connected to Z(B) via 

F(T,fx) = -T\n(Z(B)) B . (6) 

It is clear that ([6]) contains all possible interactions, perturbative and non- 
perturbative (NP) between quarks and gluons, and in particular creation or 
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dissociation of bound states. It is impossible to take into account all possible 
subsystems and interaction between them, and it is imperative to choose the 
strategy of approximations for the quark-gluon plasma as a deconfined state 
of quarks and gluons. 

We assume that the whole system of N g gluons and N q , Nq quarks and 
antiquarks for T > T c stays gauge- invariant, as it was for T < T c , however 
in case of deconfinement and neglecting in the first approximation all per- 
turbative and NP interaction, any white system will have the same energy 
depending only on the number and type of consistuents. In this case Z fac- 
torizes into a product of one-gluon or one-quark (antiquark) contributions 
and we calculate the corresponding thermodynamic potential. 

This first step is called the Single Line Approximation (SLA) and in the 
next two subsections we calculate the known results for free gluon and quark 
gas in our path-integral formalism following [10J. 

However already in SLA there exists a strong interaction of a gluon (or 
q, q) with the NP vacuum fields. It consists of colorelectric (CE) and col- 
ormagnetic (CM) parts, as shown in section 3. The CE part in the decon- 
finement case creates the individual NP self-energy contribution for every 
gluon or quark, and this is a factor of the corresponding Polyakov loop dis- 
cussed in section 4. Note, that the quark (gluon) Polyakov loop factor is 
computed from the gauge invariant qq{gg) Wilson loop, which for the NP 
Df interaction separates into individual quark (gluon) contributions equal 
to the modulus of the corresponding Polyakov loop. Therefore our Polyakov 
loop factors are expressed through color singlet 9heavy-quark) free energy 
F s (oo,T) and the popular Z(N C ) symmetry is irrelevant for them. 

With the CM part the situation is more subtle. Strong CM fields in the 
deconfined QCD vacuum introduced in [TU], [T7| and measured on lattice 
[T51 [H] create bound states of white combinations of quarks and gluons, and 
this is formally beyond SLA. To take these bound states into account one 
needs to use source terms, introduced in Appendix 1, and do calculations as 
shown in Appendix 2. 

Therefore the CM fields are taken into account as two-body correlations, 
which are discussed in section 5 and Appendix 2. 

The SLA result for the pressure with the CE is given in section 6. The 
effects of possible qq, 3q,gg,3g, etc bound states which might be important 
nearby T > T c , as it is discussed in [23], are deferred to future publications. 
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2.2 The lowest order gluon contribution 

To lowest order in ga^ (keeping all dependence on gB^ explicit) we have 

Z, = e~ F ^' T = N'(exp(-F (B)/T)) Bj (7) 
where using Eq. (jSJ) F (B) can be written as 

^F$ l {B) = ^lndettr 1 -lndet(-D 2 (5)) = 

= 4;r f(1 )V t r{ (! )V«}. (8) 

y 1 Jo s Jo s J 

In Eq. (jSJ) Sp implies summing over all variables (Lorentz and color 
indices and coordinates) and £ (t) is a regularization factor. Graphically, the 
first term on the r.h.s. of Eq. (jSJ) is a gluon loop in the background field, 
while the second term is a ghost loop0. 

Let us turn now to the averaging procedure in Eq. ([7]). With the notation 
/ = —Fq(B)/T, we can exploit in Eq. j7j) the cluster expansion [26] 

(ex P /) B =exp(£«r»^) 

= eM(f)B + l[(f 2 }B-(f)l}+0(f)}. (9) 

Below in this paper we consider mostly the leading term in the cluster ex- 
pansion, namely (/)#, leaving the quadratic and higher terms for the section 
7. 

To get a closer look at (<p)b we first should discuss the thermal propa- 
gators of the gluon and ghost in the background field. We start with the 
thermal ghost propagator and write the Fock-Feynman-Schwinger Represen- 
tation (FFSR) for it [10] 



_ roo _ roc _ 

i-D 2 )~i = (x\ I dse sD2 ^\y) = jf d8{Dz)%,e- K *(x, y ). (10) 

Here $ is the parallel transporter in the adjoint representation along the 
trajectory of the ghost: 

$(x,y) = Pexp(ig J B^dz^) (11) 



1 In what follows we shall call the contribution of gluon (ghost) Green's functions 
G(D~ 2 {B)) the gluon(ghost) line contribution to distinguish it from the loop expansion 
in standard perturbative analysis. 
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and (Dz)™ is a path integration with boundary conditions imbedded (de- 
noted by the subscript (xy)) and with all possible windings in the Euclidean 
temporal direction (denoted by the superscript to). We can write it explicitly 

as 



(Dz) 



■HI 



ff -» m =i (4vr^) 2 
d A p 



E 



cxp 



n=0,±,.. 



N 



\m=l / 



(12) 



Here, ((k) = z(k) — zik — 1), Ne = s. We can readily verify that in the free 
case, = 0, Eq. (fTUl) reduces to the well-known form of the free propagator 



ds exp 



1 



4s 



(27T) 



x exp 

E 



f° exp I" 



-p 2 s 



with 



d((m) = 



— ip(x — z) — ip^n(5 
dC(m) 

(47T£) 2 ' 



ds 



d A p 



(13) 



6(p 4 /3 - 2nk) 



Using the Poisson summation formula 

— J2 exp(ip 4 n/?) = 

An n=0,±l,±2... fe=0,±l,... 

we finally obtain the standard form 

Td?p exp[-ipi(x - y)i - i2rckT(x 4: - 2/4)] 



-d 2 )- 1 



E 

fc=0,±l. 



(2*) 



pf + (2^kTf 



(14) 



(15) 



Note that, as expected, the propagators ffTUj) and (fl5j) correspond to a sum 
of ghost paths with all possible windings around the torus. The momentum 
integration in Eq. flT2J) asserts that the sum of all infinitesimal "walks" 
C(m) should be equal to the distance [x — y) modulo iV windings in the 
compactified fourth coordinate. For the gluon propagator in the background 
gauge we obtain similarly to Eq. ([TP]) 



G 



ds(Dz): y e- K $ F (x,y), 



(16) 
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where 

rt 



$ F (x, y) = P F Pex P (2ig jf F(z(r))dr^j exp [ig J* B^dz^ . (17) 

The operators PpP are used to order insertions of F on the trajectory of the 
gluon. 

Now we come back to the first term in Eq. (jHJ), (<p)b, which can be 
representated with the help of Eqs. ffTOl and ffT6]) as 



(Ff(B))B = -Tj ^{s)d*x{DzY xx e- K 



~tr($ F (x,x)) B - (tr$(x,x)) B 



where tr implies summation over Lorentz and color indices. One can easily 
show [10] that Eq. ffTBl) yields for = the usual result for the free gluon 
gas: 

Ff(B = 0) = -Tf(B = 0) = -(N 2 - l)vX^. (19) 



45 

2.3 The lowest order quark contribution 

Integrating over the quark fields in Eq. (j3J) leads to the following additional 
factor in Eq. (J5]) 

det(m g + D(B' + a)) = [det(mj - D 2 (B' + a))] 1/2 , (20) 

where B'^ = B^ — ^-5^, fi q = \x is the quark chemical potential, and m q is 
the current quark mass (pole mass when next terms in a s are considered). 
In the lowest approximation, we may omit a M in Eq. (|20|) . As a result we get 
a contribution from the quark fields to the thermodynamic potential 

IfS{B>) = ilndetK -D\B>)) = -\ S p f £(s)^ e~ m (21) 
where Sp has the same meaning as in Eq. (jSj) and 

D 2 = (D^,) 2 = D 2 (B') - gF, v a, v = D 2 — gaF; 



1 



a 



fJ,U 



1/ 



{lulu-lulu)- (22) 



Our aim now is to exploit the FSR to represent Eq. (|21| in a form of the 
path integral, as was done for gluons in Eq. ffTUj) . The equivalent form for 
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quarks must implement the antisymmetric boundary conditions pertinent to 
fermions. We find [TU] 



1 I Jo s 

where tr implies summation over spin and color indices, 

W a (C n ) = P F P A exp (ig f B'^dzA exp g f io uv F av ) dr, 

\ J C n / JO 



(23) 



and 



(Dz] 



xy 



X 



£ (-1; 



d A p 



exp 



N 



■ / N 

ip £ c( 

\m=l 



m) — (x — y) — 



(24) 



n=0,±l,±2,... 

(2tt)^ 

It can readily be checked that in the case — 0, fi = the well known 





expression for the free quark gas is recovered, i.e. for m q 

7ir 2 

' N c V 3 T A -n f , 



Fq (free quark) 



180 



(25) 



where rif is the number of flavors. The derivation of Eq. (1251) starting from 
the path-integral form ([231) is done similarly to the gluon case given in the 
Appendix of the last reference in Ref. [10] . 

The path C n in Eq. (|23|) corresponds to n windings in the fourth direction. 
Above the deconfinement transition temperature T c one sees in Eq. ( 1231) the 
appearance of the factor 



n 



fund 



Pexp 



19 



In the next section it is shown that it will generate a gauge-invariant quantity 

Lfund = -^-(tfQfund) ■ (26) 

where Lf un d is the Polyakov loop in the fundamental representation. 
For a nonzero fi and rif massless flavors with B u = one has 



(F) 



7T- 



-™/£ 



\n+l 



-cosh 



n=l 



n' 



(27) 



In what follows, we consider arbitrary fields B n (x) and express (F) through 
field correlators. 
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3 Nonperturbative dynamics of single quark 
and gluon lines 

As was shown in the previous section, Eq. (18), the pressure for gluons can 
be written as (the T-independent term with n = is subtracted, and the 
relation P g iV 3 = —(F (B)) B is used), 

P 9l = (N 2 c -l) r^Y,G M (s) (28) 

and for quark one can write using (123]) 

P q = 2N C / — e~ m > ^(-l) n+1 [S (n) (s) + 5 { " n) (s)]. (29) 

^° S n=l 

Here and are the proper-time kernels of gluon and quark Green's 
functions at nonzero T, which are defined as follows 

GM{s) = J (Dz)Ze- K tr a (W s (C n )) (30) 

S {n) (s) = J WzT on e- K tr f {W a {C n )) (31) 

where 

ir f W a (C n ) = ^tr f W a (C n ) (32) 

and 

tr a (W s (C n )) = J^TY^FiCn) - *{C n )) (33) 

In what follows we shall neglect the first exponent in (17) both in W-£ 
and in W a , since this gives spin-dependent interaction, which can be treated 
as a correction (we relegate this treatment to the Appendix 2). Therefore 
the calculation of G^ n '(s) and S^ n \s) needs first the evaluation of the Wilson 
loop (W(C n )), 

Our purpose is to calculate G^^s) , S^(n) using for (W(C n )) the field 
correlators. 

To proceed one should look more carefully into the topology of the Wilson 
loop Wz >CT (C n ), which is to be a closed loop in Ad. At this point one should 
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again emphasize that the basic states \k > with the quark-gluon numbers 
Ng,N( k \Nq entering into the partition function Z are gauge invariant, 
Z = Y,k(k\ e ~ H ^ T \k) an d one write for \k) a generic decomposition of the 
type, \k) = \(gg)(gg)(ggg)(qq)---) where particles in parentheses form white 
combinations. 

Correspondingly (and neglecting interaction between white subsystems) 
for each white subsystem one has contribution proportional to the product, 
e.g. for the (gg) or (qq) system 

Z (la) = J JdT.dT^trW^C^,^)}, 

with 

J dT l = J d Si e- Ki {Dz®). 

Here (trW^ a (C^C^)} is the closed Wilson loop made of paths (C^C^) 
of two gluons (or qq) and parallel transporters (Schwinger lines) in the initial 
and final states, necessary to make them gauge invariant. 

The important result of this paper, shown below, is that colorelectric 
correlator Df yields factorized contribution, i.e. Z^ 2 ) —> ZiZ 2 , where each 
of the factors Z t contains only a part of the common loop in the form of the 
Polyakov loop factor (with the singlet free energy in the exponent). 

The contribution of colormagnetic fields does not factorize, however, as 
shown in Appendix 2, this contribution can be considered as a correction 
therefore in this section it will be neglected, leaving the topic to section 5. 

As we shall see, the 4d path integral will be decomposed into a product, 
D A z — > D 3 zDz^, and the interaction kernel also factorizes in the Gaussian 
approximation (W(C n )) = (Wi)(W 3 ), which leads finally to the relation (42). 
To proceed we shall use as in [23] the Field Correlator Method, where for any 
closed loop C n one can apply the nonabelian Stokes theorem and Gaussian 
approximation and write (see [6j [7] for discussion and refs) 

fr f (W(CPC®)) = ^(Pexpig f B,dz») = 

1 V c J C n 

= ^exp (-£ jf J s da^(u)da Xa (v)(F^(u)F X(T (v)) \ . (34) 

Here one should specify the surface S n with the surface elements da^ v {u) in 
the integral on the r.h.s. of (1541) . For this analysis it is important to consider 
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different terms in the sum over //, v in (I34p separately. For the field correlator 
in (34) one can use the decomposition (second ref. in [5]) which for nonzero 
T should be written separately for colorelectric field Ei and colormagnetic 
field Hi [T0l[T7]. With the definition H k = \sijkFij one has 



g 2 (ir f [Ei(x)${x, y)E k {y)${y, x)}) 



8 



ik 



D E + Df 



2 dDf 



+ UiU k 



dDf 
du 2 



g 2 (ir f [H i (x)^(x,y)H k (y)<!>(y,x)}) 



8 



ik 



D H + Df + u 



H 



H 



du 2 



g 2 (ir f [Ei(x)^(x, y)H k (y)$(y, x))) 



dDj 
du 1 

dD» E 



■ UiUk- 

1 

2 OU4 



(35) 



We take first the term with v = o = 4 (note that by definition \x < v and 
A < <7, see jSld). 



J E = - 



Da, k 4{u, v)da i4 (u)da k4 ,(v) 



(36) 



where we have defined D^xa = g 2 {ir /i 7 ^ («)$(«, v)F\ a (v)$(y, u)) . At this 
point one takes into account that the surface S n is inside the winding Wilson 
loop for the gauge-invariant qq system, as it is done in [21] and Appendix 
2. The contribution of Di leads to the sum of terms for q and q separately 
of the form (I3"8l) . which fact justifies the use of singlet Vi(oo) for single line 
quark (or gluon) term in P g (or P g i). For a single quark one obtains in this 
way 

1 rn/3 rn/3 roo > 

Jn=-J du, / dv 4 / ^f(VK^W) = 

2 Jo Jo Jo v 



J"tdZD?(y/v* + e 



(37) 



Note that for n = 1 one recovers the expression for the Polyakov loop ob- 
tained in [23], namely 



'fwd = exp(-— Vi(oo)), —Vx(oo) = Jf 



(38) 
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At this point an important simplification occurs. Namely, originally the path 
C n is a complicated path z^{j) in 4d, however the integral (37) over is 
always from some point z^ir) on C n to infinity and does not depend on the 
exact form of C n , and is the same as for the straight-line Polyakov loop. This 
is true, however, only for the -Dj4.fc4 and not for D^im- 

Another important observation is that correlators D^^iu^v) are not 
periodic function of (m 4 — t> 4 ), in contrast to fields F^ u {u) and F\ a {v). This 
will be true also for path integrals from any point x to arbitrary point y, 
(temperature Green's function) and is a consequence of the vacuum average 
of parallel transporter $(it, v) present in D^^iu.v). 

We turn now to other contributions in (I34p . 

1 f 

Jn H = 2 / D nM( u i v )dviA{u)do k i(v) . (39) 

This term is treated in Appendix 2, and is shown to be T independent, 
we shall neglect it in what follows. 

We finally turn to the spatial term, which will play in what follows a spe- 
cial role. Here the main contribution comes from the colormagnetic correlator 
D H (u — v), which provides the area law for the (closed) spatial projection A% 
of the surface S n . Correspondingly we shall denote this term as (for A 3 3> A 2 , 
where A is the gluon correlation length, D (x) ~ exp(— x/X)) 

(W 3 (C n )) = exp(-a s A 3 ) (40) 

where 

°s = \j D H (x)d 2 x (41) 

and A 3 is the minimal area of the spacial projection of the surface S n . As in 
the Jjf case, one should start with white states of qq, gg or 3q. As it is shown 
in Appendix 2, the colormagnetic vacuum (D ) acts only in the L ^ states 
of these systems, and the resulting contribution of D H does not separate 
into single-line terms, but rather acts pairwise (triple-wise for 3q). Therefore 
one can account for color magnetic interaction in the higher (2-line or 3-line) 
terms. One can take it approximately averaging single line in the gas of 
quarks and gluons, and denoting the corresponding term as Magnetic Loop 
Factor (MLF) Gs(s) and 5*3 (s) for gluons and quarks respectively. In what 
follows we shall neglect colormagnetic interaction in the first approximation. 
As a result the 4d dynamics in (1301) . fl3~ll separates into 3d and Id, and one 
can write 
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G M( a ) = G { r\s)G 3 (s) = J (Dz 4 ): n e- K ^G 3 (s) (42) 
and similarly S^(s) = S^ 1 \s)S3(s), with the free quark (gluon) factors 

Si°\s) = Gf\s) = J (Z^)ooe-^ = ^_L_ (43) 

One should notice at this point, that G^^s) and S^ n \s), Eqs. (30) and 
(31) differ in spin-factors (17) and (23), and in color group representation. 
Dropping spin-factors in the lowest approximation, one obtains from G^ 
replacing adjoint quantities (like a s (adj)) by fundamental ones (a s (f)). At 
this point one can use Casimir scaling [231 [21] to write a s (adj) = ^a s (fund). 

To compute G±(s) one can notice, that does not depend on z 4 and can 
be taken out of the integral, while the integral over z^ can be taken exactly, 
namely splitting the proper-time interval Ne = s, N — > oo one can write 



d l± eV ^- v \ s = _J_ e -^ (44) 

2 71" 2 \/ ITS 



As a result one has for G^ n Xs] 



G <»)( S ) = * c -=£-#Gf (a) (45) 



Z\/ TVS 4 



In the next section we consider the Polyakov loop factor Jjf 



4 Properties of Polyakov loops 

As it was discussed in the previous section, Eq. (128jl . the lowest order contri- 
bution to the free energy of the quark and gluon contains the path integral 
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S^(s) and G^(s), Eqs. (13"Uj) and (I3"T|) . which is expressed through the 
Polyakov loop with n windings, namely (cf Eq. (j45p ). 



G^\s) = ^=e-^Gf\s)L^ (46) 



where 



Lg. ee exp(- J* ), J,f = ^ (47) 



and Jjf is given in ([3 



n/3 [ n P , f v 



In what follows we shall discuss the following properties of Polyakov loops: 

a) Dependence of on n. (in what follows L stands for Lj, i = fund, 

adj.) 

b) Renormalization of . 

c) Calculation of L a dj and Lf un d through the field correlators, 
a) We start with n — 1 and write as in (|38l) 



= exp (-— ^(oo)J , —^(00) = ^ (49) 

where V\ (r) was introduced in |23j as a potential between heavy quark Q and 
antiquark Q at a distance r, obtained from the correlator of the Polyakov 
loops. As it is seen in (T4"51) Vi(oo) (as well as Vi(r)) is calculated directly 
through the correlator Df (for T > T c ). The latter has perturbative and 
nonperturbative contributions which we write as (we neglect posssible inter- 
ference terms, which contribute near r = 



Df(x) = Df pert (x) + Df nonp {x) (50) 

and it is clear that Df ert contributes to the color-Coulomb law between 
Q and Q and to the perturbative self-energy correction, which should be 
properly subtracted and renormalized, as discussed below in the point b). 
The NP ingredient to L is given by Df nonp and this is responsible for the 
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NP thermodynamics of the qg plasma, which is the leading contribution 
considered in this paper. 

The form of the Df nonp has been found previously at T = analytically 
[TB] [23| [18] and on the lattice [131 El an d it was shown that Df nonp can be 
connected to the gluelump Green's function, the gluelump mass Mi, found 
in [2E], defining the range of Df nonp , 

Df nonp {x) « const exp(-Mi|x|), |x| > l/M l (51) 

In [T3J E] -Df and the value of M\ was computed on the lattice, and 
in [23] this calculation was compared to the lattice data [29] for Polyakov 
loop correlator showing a reasonable agreement. Therefore one can use as 
an acceptable the value of M\ « 0.7 4- 1 GeV, which is comparable to the 
inverse of gluon correlation length AatT = 0; 1/A ~ 1 GeV. 

The asymptotic behaviour (51) is also kept at T > T c , as shown by lattice 
data P31CEI]. 

At this point one should use the behaviour (|5"T1) in (jUJ). It is important 
to stress, that Df is not periodic in v = £4 — 2/4, we remind that 

Df ~ {trE l {x^{x,y)E l ^y)) (52) 

in contrast to E^x), which satisfies PBC E^x^ + n(3, x) = _Ej(x 4 ,x), and this 
is a consequence of the presence of the factor $(a;, y) which is not periodic in 
X4 — 7/4. Hence Df(i/,£) is decreasing for large \v\ exponentially as in flBT]) . 
One can define in the integral (HHl) two possible regimes: 1) /3Mi ^> 1, and 
2) ^Mi < 1. 

In the first case one can always neglect in fHHj) all dependence on n except 
for the first factor, writing 

Jn « y r dp r ^°^ 2 + v2) - nj °- (53) 

This results in a simple relation 

4 n) = (Li) n , i = adj, fund. (54) 

In the second case Df extends over the region of several periods n, and 
for such n, that n(3Mi < 1, behaves quadratically in n. Since Df is 
positive, this means that (L n ) decreases at large T,T> M 1; faster than is 
given by the law ([53]) . Hence the general conclusion is that for T Mi the 
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law (1511) holds, while for T 3> Mi (151)) is replaced by a stronger dependence 
of the type »s (L) n for the first n terms satisfying n < n* = -j^. 

b) We now turn to the point b) and discuss renormalization of Polyakov 
loop, which was introduced on the lattice in [29] and was considered analyt- 
ically in [23J. 

In [23] it was found that the Polyakov loop expressed through the corre- 
lator Df, can be written as 

Lfund = exp( — ) (55) 

r (np) , 



where Vi (oo) is expressed through the NP (and regular at x = part 
T ^inon P ( x )) > while Vi ert \oo) is expressed through Df t (" 1 " r ' ~ n ■ 1 >,, i f - " 



the lattice-type minimal distance a in the lowest order in g is equal to 

V (pert) (oo) ^^a 1 (l_ Tlna \ (56) 

n \a J 

Neglecting the next order and possible perturbative-nonperturbative contri- 
butions, one can define the renormalized Polyakov loop as 

^ exp (_&>) (57) 

This definition can be compared to the renormalization of color singlet 
free energy F s (r) used in [29|, where it was noticed that F s (r) does not 
depend on temperature (at least for T < 2T C ) for distances r < 0.15 fm. 
Thus normalization of the data at short distances allows to have a smooth 
limit to the T = case. 

Since V np (r) ~ r 2 [23||5] at small distances, F s (r) actually coincides there 
with Vx P (r), and the divergent part of Vf ert \r) is in Vi P (oo). Hence 
one can deduce, that normalization to the small distance part of V s (r) at 
T = 0, with 

F r s en (r) S V x (r) - V} pert) (oo) = v} ren) (r) (58) 

is a justifiable renormalization procedure equivalent to dropping the divergent 
part, Vi Pert \oo) in the definition of Vi ren \ In what follows we shall define 
L ren as in (JS@. 

One might ask how at all the (unphysical) - divergence can occur in the 
expression for the free energy (pressure), where physical quark and gluon 
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Green's function enter. Indeed, we have used the path-integral form (I30|) . 
(J3TJ) , which is exact and corresponds to the standard perturbation theory at 
small distances. Therefore self-energy parts due to gluon exchange have a 
normal logarithmic divergence which is eliminated by a standard renormal- 
ization procedure, setting the renormalized self-energy equal to the pole mass 
for quarks or zero for gluons. Therefore the divergent part is effectively elim- 
inated from the Green's function when general form as in (T2SJ) is used. Note, 
that the factorized form (142|) already used the fact that Jjf is a constant, 
not depending on z&, z, which implies that the self-energy divergence is elim- 
inated and only nonperturbative contribution to L is kept. Therefore Jf in 
( 145]) is already Vi np \oo)/2T and corresponds to the renormalized Polyakov 
loop flSTj) . 

At this point one should be careful to distinguish between the (theoretical) 
definition of the Polyakov loop (57) and lattice values L\ at , since the latter, 
as discussed in [23] , contain all possible excitations of the heavy-light system, 
not present in the light quark or gluon loop. For more discussion see last 
section of this paper. 

c) Here we discuss calculation of L f und and L ad j through D and D 1 both 
below and above T c As found in [23], L ad j is 

L adj = exp {- (~V D (r*,T) + ^(oo.T)) /t} (59) 
where Vd is given in Appendix A of [23], T <T C 

V D (r, T) = 2 f dv{\ - uT) f (r - S)d£D E (J? + z/ 2 ) (60) 
Jo Jo v 

and r* is an average distance between the adjoint line (static adjoint charge) 

and a gluon of vacuum, forming together a bound state (a gluelump). In 

a rough approximation one can replace |yb(r*,T) ~ M g i p , with M 9 ; p -the 

gluelump mass of the order of 1.1 4- 1.4 GeV [28]. In the similar way one has 

for a quark Polyakov loop 

L fund = exp{-(V D (rt,T) + i^(oo,T)))/T}, L adj = {L fund f' A (61) 

where r* corresponds to the average size of the heavy-light meson, r\ ~ 
1.2 4 1.4 fm, and V D (r*,T) w e HL w 0.6 4 0.7 GeV. 

Therefore both L a( § and Lf und for rij > are nonzero also below T c , 
however small in this region. Above T c one has a Casimir scaling relation, 
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derived first in [23J as a consequence of the fact that Gaussian correlators 
D, Di proportional to the charge squared. 

Recently the Casimir scaling as in (61) was accurately checked and con- 
firmed on the lattice [24] for five lowest representations of SU(3). This fact 
shows that also for Polyakov loops as well as for Wilson loops [30] the lowest 
(quadratic) field correlators D, D\ are dominant with accuracy of few per- 
cent. Hence the QCD vacuum is Gaussian dominant at all temperatures, 
(unless a finetuned cancellation of higher correlators takes place [30J). 



5 Nonperturbative Equation of State (Single 
Line Approximation) 

We are now in position to summarize our derivations of the nonperturbative 
free energies for quarks and gluons in the Single Line Approximation (SLA) 
as follows: for gluons from (28), (45), 

sla _ CjV| - 1) r ds n®( B \ ST r W 



(iv c 2 - 1) r°°ds^ _je 

1 2_^ e J ™ (62) 



16tt 2 Jo s 3 n ^ Q 



where J® is given in (47), (48), (53). 

For quarks of rif flavours from (29), (45), 



v 71 " Jo s ' n=1 1 

N c f 00 ds 2 ™ ,, JE an . . 

= nf 4^l ~s* 6 9 E(-!) e ^ e n cosh^, (63) 



n=l 

The Polyakov loop exponents, , are expressed through the function Df in 
(37) and can be either computed analytically, as in [23], [18], or taken from 
the lattice data (see e.g. [3~T]). 

As a first check we consider the limit of free gluons and quarks, i.e. J„ — > 

1. 
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Integrating over ds in (1621) . (1631) one obtains for the free gas with m q = 0, 
H = 



P s1V 0) = _ 0, = Z^„ (64, 



where we have used the standard sum values 

Vl = - V ( ~ 1)n+1 - 77r * f65) 
^n 4 90' ^ 720' 

Another simple approximation consists of neglecting all terms in the sum 
over n, except the term with n — 1. 

The accuracy of this approximation can be estimated from the sum (165j) , 
where all terms with n > 1 contribute less than 8%, and for nonzero Jjf 
(for T < 1.5T C , Lf un d is less than 0.7 for n/ = 2 [32J) higher terms in n are 
suppressed even more. 

Therefore with accuracy better than 10% (for ^ < 1) one can suggest the 
following expressions for -Pj LA , P^ LA 

P 9 f A = ^f^T 4 (66) 

P q SLA = ^ t tl) ^(Lf U n d r^\T)cosh^ « ^^n/L /Mnd T 4 ^W (T)cosh 

(67) 

where we have defined 

A representation similar to (1661) was suggested a decade ago in [12], for 
gluons where instead of as in (1661) . the Stefan-Boltzmann factor |g was 
used. With the lattice data for L ad j known at that time, the resulting agree- 
ment of Pgi LA was quite satisfactory. 

In the next paper [33] we exploit new accurate data for the renormalized 
Polyakov loops, fundamental with rif = 2 [31] and rif = [29J and recent 
data for the adjoint loop [24J, showing a reasonable agreement with fl66|) . 

§3). 

Summing up the series ( |67l) for m q = one easily obtains another form 
for P g SLA 
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^ qrA/ . n/T 4 [ r 00 z A dz r°° z s dz ] , , 

P q SLA (m q = 0) = J—\ + (69) 

and and for m q ^ (see [33] for details) one has instead 



(70) 

where z/ = ^r. 

Another useful quantity to compare with lattice data is the internal energy 
density e, 



and from (1661) . (IBTj) we obtain (for = 0) 



4 = ~2 T ( N c ~ 1)^(^0 (72) 



4 LA = -^n/r 2 — (I*L /un ^Ci)(r)) (73) 



4iV c 2 d 
~i^~ nf dT 

and finally for the "nonideality" of the quark-gluon plasma one obtains 

e SLA_ S pSLA 2 {Nl-l) m d (T 

+ ^ T ^ (L/W ^ 1)(T)) - (74) 

It is interesting to study the limit of large T for ' T4 given in Eqs. (J66|) . 
( IBTj) . and compare it to the Stefan-Boltzmann value ^r(T = oo). For that 
purpose one needs the large T limit of Lj, i = gl, q- 

The situation is more transparent with the large limit of ip^(T). 

From (81) one can conclude that ip)p(T) tends asymptotically to 1 with 

- 2 

the (small) correction proportional to j§ . As to L iy from (47), (48), (49) 
it is clear that the problem with Lj boils down to the large T behaviour of 
Vi(oo, T)/T. In this paper we take L, as given by lattice [271 [291 EI] and 
analytic calculations, [151 (HI 121] from which Lj grows with T and reaches 
the value of 1 around T = 2T C , and may become larger than 1 beyond 
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3T C . This is connected to the negative value of the lattice free energy of 
static quark in this range of temperature [32]. The actual limit of Lj at 
infinite T is connected to the fact, how the magnetic confinement, viz. cr s (T) 
enters in Li, and this will be discussed in one of the next papers. In this 
way one expects that all corrections to — 1 should be proportional to 

some power of ^ a ^ T ^ ~ O ( lnT / A ), since a s (T) is the only nonperturbative 
dimensionful quantity which defines the dynamics of high temperature QCD 
- colormagnetic interaction [TTJ [19] . 



6 Beyond Single Line Approximation 

In previous sections one-particle contributions were considered, i.e. SLA was 
used, however it was stressed, that some interactions like the colormagnetic 
one, acts in white ensembles like qq, gg etc. These latter contributions gen- 
erated by the source terms listed in Appendix 1 were neglected except for 
contributions of Df , which acts in qq(gg) but produced effectively single-line 
term Vi(oo). In addition Df as well as D H , D 1 ^ produce qq, gg, 3q etc bound 
systems [23j [18j [M] which may affect strongly dynamics at T > T c , as will 
be discussed below. 

In general one can write the contribution of the pair (triple) interaction, 
starting from the pair (triple) Green's function, e.g. for the qq or gg system 
as 

Gv2 = J dV.dV^trW^C^, Cf )) (75) 

with dTt = f™ d Si {Dz)l m e- K *. 

Here the closed Wilson loop is formed in the same way, as discussed at 
the beginning of section 3. From the (qq) or (gg) Green's function one can 
proceed to the derivation of the two-body Hamiltonian H q q or H gg in the 
same way, as it was done for T = in [35j, 

G u (x 1 x 2 | 2/i2/ 2 ) = (xiX 2 \e~ Hl2/T \y 1 y 2 ) (76) 

Our purpose in this section is to discuss the properties of the resulting 
Hgg, Hgg for T > T c , and possible consequences for the quark-gluon thermo- 
dynamics. Writing Hi = Ki + Vi, i = (qq), (gg), (qg), one can consider 
two distinct dynamics, the Colorelectric (CE) one due Df , and colormagnetic 
(CM) due to D H , D*. 
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In the CE case the interaction can be written as (see [2"3"l [IB"] for more 
discussion) 

V q f E \r,T) = Vl (r,T) = V^T) - ^(oo,T) (77) 

and [23] 

Vi(r,T) = / VT (1 - vT)dp [* td£D*{ifi* + ?) (78) 
jo Jo 

In (ITT)) we have subtracted Vi(oo,T), which is already accounted for in 
SLA in the form of Polyakov loops, see Eq. (55). 

Bound states of qq,gg, and gq systems have been studied in the field 
correlator formalism in [231 ITB"] and on the lattice (see [29] and refs. in [23J , 

M)- 

Characteristic feature of v i(r, T) is that it is short-ranged (r e // ~ 0.3 fm) 
and can support bound S-states in (cc)i, (gg)i(cg) 3 , (gg)s of weak binding, 
|e| < 0.14 GeV for T < 1.5T C . 

On the lattice in addition to 5-states of (cc)i, (66) i and light (qq) also 
lowest P state is claimed to exist (see [36] for a review). 

It is clear that weakly bound states should fast dissociate in the dense qg 
plasma and therefore hardly produce any significant effect on the production 
rate, however can be important for the kinetic coefficients. 

We now turn to the CM case. Here dynamics of qq system was carefully 
studied in |34j using the formalism of [35]. The resulting interaction has 
several interesting features. 

1. Interaction of quark and gluon systems with the CM vacuum occurs 
only in the states with the orbital momentum L > 0. 

2. Centrifugal barrier is effectively killed for r 2 3> L ^ +1 ^ and no spin- 
independent interaction is found in this region. 

3. Long-range spin-orbit forces are predominant at large r, 

V ao {r) = — LS(i s / filjj-r, where y? e ^ is the effective mass squared to 
be determined from the extremum of the relativistic Hamiltonian (we 
disregard in this section the large r damping of V so (r) due to string 
inertia in \i. 

The latter can be written as (the qq system with equal masses is treated 
here for simplicity) 

2 2 

H u = V -i±^ + n + v<g>(r) + V so (r) (79) 
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where the sum of potentials can be effectively rewritten as (see [M] for exact 
expressions) for J = L + S 

L(L + 1) o~cL 

V C m{t) = ^ - '- . 80 

[fi + o s r)r l 2n(n + or)r 

Note that Vcuij') = for L = 0, as was noted before. For large r the first 
term decreases as O (jg) . 

Numerical analysis shows that the lowest bound states appear in the 
situation when r < fi/a, when Vcm{^) can be rewritten as 



V M~ L ( L + 1 ) a * L L(L + l) a eff 
VcM{r) = 5 = 5 (slj 



with« e// = ^|. 

Note that is to be found from the minimum of the Hamiltonian eigen- 
values, as it is always done in the einbein Hamiltonian formalism [35]. More- 
over VcAi(r) and especially V so {r) is not treated as a correction: in contrast 
to the standard expansion, the spin- dependent interaction in the Field- 
Correlator Method is obtained in the framework of the Gaussian (quadratic) 
approximation, which is proved to be accurate, see [30], [21] for more details, 
and [37] for the derivation of spin- dependent interaction. 

For the fixed fi, one obtains the bound state mass from Hu as 

^M=-+»- 2{n ;j> +ir (82) 

The physical mass of CM bound (qq) is obtained after an extremum in \x 
is found for some [i = fiQ. As was found in [34J, this extremum exists for 
m > 0.25 GeV, i.e. for (bb), (cc) and possibly (ss) systems with the binding 
energies 0.007, 0.19 and 90 MeV for a s = 0.2 GeV 2 . 

There is no extremum for m < 0.2 GeV, which means that relativistic 
system of light qq is not stable (no lower bound for the mass). The situation 
is the same as for the famous Z > 137 problem in QED, where relativistic 
treatment of bound states for a e ff > 1 leads to inconsistencies (see [3S] 
for reviews). Physically the problem should be treated taking into account 
pair creation, i.e. in the one-body (heavy-light) Dirac formalism, as in [38J, 
considering also lower continuum of negative energy states. 

In the framework of the relativistic (einbein) Hamiltonians [35J this be 
studied using the matrix Hamiltonian, suggested in 
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Thus light quarks in the (qq) system or together with heavy antiquarks 
are unstable in the CM vacuum with respect to multiple light pair creation. 
A similar, and even more pronounced situation should occur for gluonic sys- 
tems, like (gg) or a gluon circling around a heavy quark. Here Hamiltonian 
is the same as in f!79|) with the replacement a s — ► |cr s and quark spin by the 
gluon spin, so that 

9 

a e ff (gluon, /i) = -« e // (quark, /x). (83) 

Again the gluon systems are unstable and should be stabilized by creation 
of multiple (gg) pairs, which means that the CM vacuum stimulates numerous 
(qq) and (gg) pair creation. 

Numerical estimates of these processes and spectra of produced quarks 
and gluons are of immediate interest for the ion-ion collisions and shall be 
treated elsewhere. 

At this point one should consider the effect of quark-gluon medium on the 
existence of CM bound pairs. It is clear, that the size of the pairs is (from 
the minimum of Vcm and a e ff ~ 1) r e // ~ ~ 0.5 fm L 3 / 2 , and in the 

medium of density n with closest neighbor distance r$ ~ n~ 1//3 , one expects 
strong screening for r e ff > r . Hence for T > 1.1T C when the energy density 
exceeds 2 GeV/fm 3 , r ~ 0.5 fm the high-L levels are already screened. This 
might explain a fast drop of entropy S^T) beyond 1.05T C found on the 
lattice [ID], as well as large value of S^T) near T c , since large -L, large- J 
bound states might bring about a large entropy. 

7 Summary and conclusions 

We have used BPTH and field correlators to find the NP contribution to the 
quark-gluon EOS. As a result the free energy (pressure) is represented as 
a double sum in powers of a™ (T) and in the number of interacting single- 
particle trajectories. In the lowest order in a™ one obtains the purely NP 
corrections to the free (Stefan-Boltzmann) result, which contains as a leading 
term the single particle trajectories interacting only with the NP vacuum 
(background) fields - Single Line Approximation discussed in chapters 1-6. 

In the next terms of this series also NP interaction between quark (gluon) 
trajectories is taken into account in chapter 6. 
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The most interesting result of the investigation above is that in the leading 
approximation the influence of the NP background can be represented by two 
factors multiplying the Stephan-Boltzmann result: one is the Polyakov loop 
for quark and another for gluon taking into account colorelectric vacuum 
fields. 

Thus the main contribution to the pressure P(T) and all other thermo- 
dynamic quantities is given in our approach by the NP impact on the free 
otherwise trajectories of quarks and gluons. Perturbative interaction is in- 
cluded in the next orders of a™, n > 1, and yields additional correlation 
between quark and gluon trajectories (line-line interaction), which together 
with the NP correlation leads to a possible formation of quark and gluon 
bound states in the region T c < T < 1.2T C and should be considered as a 
next step, discussed in chapter 6. 

As discussed in section 5 and Appendix 2, the CM interaction produces a 
rather strong potential with the dominant spin-orbit term V so ~ — ^f 2 -, which 
for quarks with momentum p ~ fj, and - ~ p yields the plasma interaction 
ratio T g = ~ |f and T q = |r g . Both T g and T g are large and presume a 
strong interacting plasma in the temperature region, where screening due to 
high pressure is not operating. At large T the Polyakov and magnetic loop 
factors tend to unity, and the pressure tends to the Stefan-Boltzmann limit, 
as it is in the standard perturbation theory. Comparison of our results, Eqs. 
( 1621 . ( 1631) and ( |66l) . ( 1671) to the lattice data is given in detail in the next 
paper [33], where it will be shown, that a reasonable agreement is achived 
already in the situation when line- line interaction is neglected, as in (166]) . 
(|67|) in the region at T > 1.2T C , while in the narrow region T c < T < 1.2T C 
one possibly needs to take into account line-line corrections. This result is in 
agreement with an earlier investigation (see the first reference of [12J), where 
also a reasonable agreement was found for the gluonic part of the spectrum. 

At this point it is interesting to compare two distinct pictures of the 
sQGP. In the first one the main dynamics of sQGP is due to interaction 
between the quarks and gluons, which should be taken into account to high 
orders, as in [25] , with possible partial resummation of the pertubartive series, 
see e.g. |41j. It is important, that vacuum fields do not enter directly into the 
picture, however can be taken into account by adjusting to the dimensionally 
reduced theory at larger T. 

From this prospective the theoretical method suggested in the present and 
next papers of this series, is approaching sQGP from the opposite direction 
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as compared to the perturbation theory. Indeed, the main contribution is 
obtained from the vacuum (NP background) and perturbative corrections are 
treated as small at all temperatures. One important result of the proposed 
method is the calculation [21] of the Debye mass mo which is expressed 
through the spatial string tension a s , mo = co\fo~ s with calculable coefficient 
Co- A good agreement with the lattice data demonstrated in [21] gives a 
strong support to the idea, that the NP vacuum is after all the main source 
of dynamics both in the confined and in the deconfined phase. The main 
outcome is the picture of almost independent quasiparticles. 

This is in agreement with the lattice measurements of the correlation 
between strangeness and baryon number or electric charge, demonstrating 
almost zero effect for T > T c |42j . 

The method suggested in the present paper can be easily extended to find 
the phase diagram in the (/x, T) plane, which is done in the next paper of this 
series. N.O.Agasian actively participated in the first stage of this work, and 
the author is indebted to him for suggestions and discussions. The author is 
grateful for useful discussions to K.G.Boreskov, S.N.Fedorov, A.B.Kaidalov, 
O.V.Kancheli, B.O.Kerbikov and M.B.Voloshin. This work is supported by 
the grants RFBR 06-02-17012, grant for scientific schools NSh-843.2006.2 
and the State Contract 02.445.11.7424. 

This work is done with financial support of the Russian Federal Agency 
of Atomic Energy. 



Appendix 1. Background Perturbation Theory for T > 

We start with standard formulas of the background field formalism [8, 9j 
generalized to the case of nonzero temperature. We assume that the gluonic 
field can be split into the background field and the quantum field a M 

A„ = B„ + a„, (Al.l) 

both satisfying periodic boundary conditions 

B^(z 4 , Zi) = B^(z 4 + np, Zi), a M (z 4 , Zi) = a M (« 4 + n/3, z*), (A1.2) 

where n is an integer and (3 = 1/T. It will be convenient in what follows 
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to use as in [S] the background Lorenz gauge 

D,(B)a, = 0, D^(B) = 8,6^ - igT b ca Bl (Al.3) 

and to split the gauge transformation as follows (in fundamental representa- 
tion) 

A^(x) = U + (x)(A fl (x) + -d,)U(x), 

9 

B^x) = U + (x)(B^(x) + -d„)U{x) (A1.4) 

a^(x) = U + (x)a II (x)U(x) 

One can see, that the form of gauge condition (1A1.3[) is invariant under gauge 
transformations (1 A 1.41) . The partition function can be written as 

Z(V,T,n) =< Z{B) > B , 

Z(B) =N J D<pexp(-J dr J d 3 xL(x,r)) (A1.5) 

where <fi denotes all set of fields a M , N is a normalization constant, 

and the sign <>b means some averaging over (nonperturbative) background 
fields Bfj,, the exact form of this averaging is not needed for our purposes. 
Furthermore, we have 

L t ot(x, t) = £ ^ + Ltf n \ a M , H/+), (A1.6) 



8=1 



where 



Li = l(F a , u m 2 ;L 2 =la«W£al, 
L 3 = e a (D\B)) ab Q b ; L 4 = -igQ a {D^ a,) ab Q b (Al.7) 
L 5 = -a{D ll (B)a ll ) 2 ]L & = L int (a 3 , a 4 ) 
L 7 = -a„D^B)F^(B); L 8 = ^+(m + /i 74 + D{B + a))* 

Here 0, are ghost fields, a- gauge-fixing constant, L 6 contains 3-gluon- 
and 4-gluon vertices, and we keep the most general background field B^, not 
satisfying classical equations, hence the the appearance of Lj. 
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An additional term in (IA1.60 . L(j( n \ ...) serves as a generating functional 
for creation of white multiquark or multigluon systems, which are singled out 
at low temperatures. One can write 

L(j (n) , a M , , ^ + ) = J dx x dx 2 {^ ] x, x u x 2 )tr(a„(x 1 )$(x 1 , x 2 )a u (x 2 )) + 

+ j { f\x, x 1 ,x 2 )tr(^ + (x 1 ) lx ^(x 1 ,x 2 )^(x 2 ))} + ... (Al.8) 

Here <&(xi, x 2 ) is the parallel transporter in adjoint representation (marked 
by tilde) made of field only 

$(x 1 ,x 2 ) = Pexp(ig [ 1 BptydZp) (A1.9) 

x 2 ) is the same but in the fundamental representaion. One can see that 
L and are gauge invariant. Differentiating in one generates white 
states of quark and gluons which enter in the total sum over n 

Z(V,T) = Y,(n\exp(-l3H)\n) (ALIO) 

n 

The inverse gluon propagator in the background gauge is 

W«l = -D 2 (B) ab • V - 2gF; v {B)r b (Al.ll) 

where 

(D x ) ca = dxSca - igT h ca B\ = d x 5 ca - gf bca B b x (A1.12) 

We consider first the case of pure gluodinamics, L s = 0, and omit the 
source term, L in (lA1.6j) . 

Integration over ghost and gluon degrees of freedom in (I A 1.51) yields 

Z{B) = A'(det W(B)); e f [det(-D^B)D^B + a))] a= _s_ x 

OO ql f 1 

x{1 + ^_M (a= _ )}exp( __ JM/ -i J)jft= D ^ (B)FMB) (A1.13) 



Appendix 2. Calculation of the Color Magnetic Interaction 
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Taking into account any interaction between gluons or quarks, and in 
particular when the colormagnetic interaction is considered, one should recur 
to gauge invariant initial and final states \l/ in the Matsubara expansion for 



P g , P q , and define ^ gg , ^ q q as follows 



■*f g \x,y) = tr(a^x)§(x,y)a u (y)) (A2.1) 

V$(a/,l/) = * gg (x + nPe i , y + n(3e±). (A2.2) 
For the double-line term in (93) one obtains 

(Z(B))%° Une = J dTdY' (trW rr >) (A2.3) 

where Wrr> is a rectangular Wilson loop, formed by parallel transporters 
Q>(x,y),<S>(x',y f ) and gluon paths (x,x') and (y,y'). 

First of all one can calculate (trWrv) using Gaussian approximation and 
cluster expansion [5], [6] . For fundamental charges (qq system) one has 

,frc Tir > / g 2 



'N, 



W rv >) = exp f-y J dir^{u)dir Xc {v){F^{u)F Xc {v))\ (A2.4) 



where we have omitted $(/u, v), &(v, u) in field correlator for brevity and 
defined 

d-n^iu) = ds^ v {u) - ia^dr (A2.5) 

and dsn U (u) is the surface element, and dr is proper time integration present 
in W a , Eq. (23). 

We start with the terms coming from the product ds^iv^dsxaiv) in 
(1A2.4j) . while all other terms have been treated in [37] and we shall write 



the resulting equation at the end of this Appendix. 

In what follows we shall derive as in [35] the effective Hamiltonian for the 
qq(gg) system, but taking into account that colorelectric string tension oe 
and the colormagnetic one <jh = &s can be different and one of them {(Je) 
may vanish. In what follows we are closely following the paper [31]. 

With the notations of (35), (36) one can write 

(ir f W$) = exp (-- J J ds^v^dsxai^D^xaiu^v)^ = 

= ex P (-J n (L>f) - J* H - J« - J n {D E )) (A2.6) 
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where J n (Df) = 2J E and J E is given in (36), while J EH is 

1 f 

J n = - E I D i4M (u,v)da i4 (u)da M (v) = 



lrr dD EH (w) 

- A J J (SikWi - S u w k ) — ^ dsu{u)ds k i{v) (A2.7) 



with w M — — v^. The integral du 4 in (1A2.7|) reduces to D EH (nf3 — V4 

n 



Df H (—V4) and is small for T <C -r^— = uMeh, moreover this integral van- 



ishes for any flat surface. 

We shall not discuss this term in what follows and now turn to the term 
J n (D E ), which can be written as 

J n (D E ) = - J J ds i4 (u)ds k4 (v)D i4M (u,v) = 

Ufa - v 4 ) 2 + fa - ^) 2 ) (A2.8) 



du i dv i duidv i D E 



1 

~ 2 

Introducing variables r, p 



r = x-y, p=[rxw], w = j3x + (1 - /3)y, (A2.9) 

where we have parametrized the Wilson surface Wrv with < t < n/T\ < 
/3 < 1, so that the surface elements are 

ds i4 (u)ds i4 (v) = Y{t)Y{t')dtdt'df3df3' (A2.10) 

while 

ds ik {u)ds ik {v) = 2p(t)p(t')dtdt'df3df3'. (A2.ll) 
One can also write using induced metrics g a b 

(u-v) 2 = (u(t,P)-v(t',f3')) 2 = g ab U b , ^ = t -t', £ 2 = /3-/3' (A2.12) 

and det g = r 2 + p 2 . In this way one finally obtains 

rn/T A r 2 

J n (D E ) = a E dt d(3^—= (A2.13) 
Jo Jo v r + P 

J n {D H ) = a H [ n/T dt ['dp j4=^ . (A2.14) 
Jo Jo v p + r z 



31 



In case when T = 0, one has D E = D , ge = &h an d one obtains 

rn/T rl i 

J n (D E ) + J n (D H ) = a dt d{3Jr 2 + p 2 (A2.15) 

Jo Jo v 

where the standard Nambu-Goto action used in [35] is Jw 2 {w') 2 — (w^w' ) 2 = 



\/r' 2 + p 2 . 

In the action for equal mass quarks, m q 
that 1 35 1 



To 



|(x 2 + y 2 )^+ / //(// 



K q + K q = 

one finally obtains the Hamiltonian for T = [32] 



mq = m, taking into account 

To 



o 



(A2.16) 



H = P r +m \p+ ^ 



2 (l + j(P-ri vdp 



where the einbein parameter v is to be found from the extremum condition. 

For the case of electric deconfinement, when oe = and <jh 7^ 0, one 
obtains instead [3i 



H 



p 2 , + m 2 
p 



+ p + ["dp 

Jo \ 2u 



v 



H V a 2 r + 

2v 2 



(A2.18) 



L(L + 1) 



2r 2 



2 



and in the case when one mass is much larger than ^an, \fo~E one obtains 
for the heavy-light system with both <jh, oe nonzero [31]. 



H, 



HL 



pi + m 2 a 

— h-+ 

2p 2 



2„2 



d/3 



erjr 

~2z7 



+ 



L(L + 1) 



2r 2 [p + J 1 d(3(3 2 u} 



where we have defined H kin = Pr ^" + §> an d 



0"! = O"// + ^(Oif - CT2 = 27]((Je ~ CT H ). 



~- Hkin + llHL 

(A2.19) 



(A2.20) 



Here rj is another einbein parameter (cf [35]). Consider first the case when 
L = 0. Taking extremum of ( 1A2.18I) . (1A2.19I) with respect to v one has 

H{L = 0) = H kin + r f l d(3[a E + (a H ~ <r E )(l ~ vf] (A2.21) 

Jo 



32 



and taking extremum in t], which yields r] = 1 one has for both Hamiltonians 

H{L = 0) = H kin + a E r. (A2.22) 

Hence for <je = 0, on ^ and L = magnetic "confining" vacuum does 
not actually confine quarks (or gluons). Hence quotation marks for magnetic 
"confinement" . 

Let us now turn to the case of L ^ 0, <je = and consider large L limit, 
when one expects as in [35] that u ^> fi , where the subscript refers to the 
extremum values of u, fi. In this case taking extremum in v in (IA2.19P and 
neglecting fi in the last term on the r.h.s. of (1A2.19I) one obtains 



?4) 



^(i + ^)V + 



3L(L + 1) 



1/2 



and for the 7] - extremum value of 7], one has an equation 



(% 4 -l)(% 2 + l) 



3L(L + 1) 



= c 



with the solution 



1 1/16 
- 3 + U{t7 + " 



c(c + 32) 



1/3 



54 



+ 



(A2.23) 



(A2.24) 



1/16 \ Vc(c + 32) 
^2 v27 + C J — 



1/3 



(A2.25) 



The resulting Hamiltonian is 



h HL = \ a%{l + V i)^ + 3L(L 2 + l) - 2 Vo a H r 



(A2.26) 



where rjQ is defined in (IA2.25I) . It is easy to check that for L = c = one 
obtains 770 = 1 and hence hn^il = 0) = 0, in agreement with (1A2.22I) for 
a E = 0. 

For small c (IA2.26I) and (IA2.25D yield 



h HL = 3L A (L \ 1 \ 3L(L + l)«alr\ 
Aa H r 6 



(A2.27) 
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In the opposite limit, c 3> 1(3L(L + 1) 3> (a H r 2 )) one should take into 
account also /i in the denominator of Hhl m (1A2. 19|) and one obtains the 
standard centrifugal term 

h HL = %£±P-. (A2.28) 
One can see, that the centrifugal barrier of light particle is essentially de- 

/ X 1 / 4 

stroyed by the magnetic vacuum for r > r , r = f 3L ^ +1 ) j while for 

r < To the barrier is essential . 

The case of equal mass particles (light-light or heavy-heavy) is done in the 
same way as above, the only difference being that instead of 3L(L+1) one has 
12L(L + 1) which should be substitued in all Eqs. flA2.23HA2.28j) . Thus the 
conclusion for the spin- independent part of the Hamiltonian, Eqs. ( ]A2.18[) . 
(IA2.19I) . is that in magnetic vacuum the centrifugal barrier is effectively 
destroyed for r > r , which fact stresses the importance of large distance 
spin-dependent interaction - to be discussed now. 

The SDI is obtained from (IA2.4I) . when one considers products of ds^ v ■ 
a p \dr ( for spin-orbit forces) and a^drapxdr' (for spin-spin and tensor forces), 
(see [37J and [43j for details of derivation). 

The resulting SDI can be represented in the Eichten-Feinberg from r,i 1 ~ 



r (diag) ( r>\ _ t °l L l _ ^2 w 1 de 2dV 1 (R) . 

sd W-Ktfi A^ 2 )[ RdR + RdR 



2/ii/i 2 R dR 12/ii/i 2 12/xi/x 2 -R 2 

where SDI potentials Vi are expressed via field correlators as in [37], see also 
definition of correlators in (35) 

1 dV l , f R dX / A 



R dR 



2 Here one should take into account that each factor dsik contains p and hence w 
and each z% for particle i yields as in (|A2.19|) the factor fk = m + J v((3)(3 2 d(3 in the 
denominator instead of /Ltj, while each brings about a factor This additional v{(3) 
term in /t, was hever taken into account before, and is inserted in [33] and below in (|A2.40p . 
dM3T|). 
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R dR 



°° , r R XdX 

dv / ~m 

-00 Jo R z 



D H (X,v) + D?{X,v) + X 



>. dD? 
dX 2 



V« = - 



00 0R 2 



V, 



roc ( f)D H \ 

J Jp f 3D H (R, v) + 3D»(R, u) + 2# 2 ^H 



(A2.31) 
(A2.32) 
(A2.33) 



1 de(R) 
R dR 



00 , r R dX 
dv I — 

-00 Jo R 



D E (X,v) + D E (X,v) + (X 2 + v 2 ) 



dDf 
dv 2 



(A2.34) 

In the deconfinement phase D E vanishes, while D H , Df , are nonzero, 
and the main part of D E , D± is perturbative, yielding the known from of 
SDI essential at small distances. One can assume (in agreement with lattice 
data [131 E]) that Df,D^ do not change with temperature, which yields 

nH _ 16a a \ 

-^1™^/ n 4 J 



K^lpert 



Ipert 3kx 



-E V 2pert( R ) 



1 



4a, _ 
3i? 3 ' R £pert 

32yra 



-^ £ 'pert( R ) - ^ 



4a, 



V4pert(-R) 



3^3' 3per 

^(R), ^ pert = 0. 



4a s 
R 3 



(A2.35) 



Since NP colormagnetic fields at T > T c act only on states with nonzero 
one can neglect V in 
(forT>T c ,/?=l/T 



L, one can neglect V^ np . Writing D E,H = D E 'f rt -\- D E ^ p ' , one has from (1A2.34D 



i2 



e' np = J (l-y/P)dyD E np ^y* + R2), 



1 z-/ 3 / rP rR XdX 

-^L v = I (l-y/(3)dyDU^Jy 2 + R2 ) + 2 1 dv{l-u/0) I -jpD H {X,v) 



(A2.36) 



Finally for V{ = V{ np one has from 
l -K = -2jf (1 - -/0*<jf § (l - 5) B »(vTO). (A2.37) 



We have corrected in ()A2.35j) . (lA2.37[) equations from [37], taking into ac- 
count the finiteness of the Euclidean time integration in the temporal Green's 
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functions at T > 0. As one can see, the presence of upper limit becomes es- 
sential only when ^ < t g , where t g is the correlation length of the vacuum, 
D H (x) ~ e" |a:|/ * 9 , for \x\ > t g , and for T = one has t g « 1 GeV" 1 . Thus 
one expects, that for T < 1 GeV~ 6T C the upper limit (3 is inessential and 
one can use the same integrals as in [37], i.e. replacing (3 by oo. 
At large R one has from (|A2.36j) . flA2.37|) 

= -a./fl, I^(i2) = ^, tf^oc (A2.38) 

where 7^ is 

/CxD rOQ 
dp \d\D H (V\ 2 + v 2 )^0(t g -a s ). (A2.39) 
-00 Jo 

Hence the main term at large R in Vs£>(A3.29) is 

V SD (R - 00) » - f ^ - ^0 ^. (A2.40) 

For equal mass particles \i\ = pi = A*, L2 = — Li = — L? S = \{(T\ + cr 2 ) 
and one obtains 

™— )-S( 1+0 (D) + fe)' (A2 ' 41) 

where A = ^ + 2 jj <W - \f v{(3) « M + ^- 

It is clear that interaction flA2.41j) can produce infinite number of bound 
states for SL > 0, i.e. for S — 1 and J = L + 1, independently of the 
dynamical mass fi. In this way (IA2.41I) embodies the suggestion made in 

One can also estimate the interaction (lA2.41j) for the free quark with the 
momentum p, when // « p and L = [r x p] = [p x p] , where p is the impact 
parameter to the neighboring antiquark. 

One has \Vsd\ ~ ~ ^fj an d the standard plasma interaction ratio Y 
for average momentum (p) ~ T 

r, = !^°U^«^; r, = 2-^ (A2.42) 
9 T fcm p 2 T 2 ' 9 4T 2 1 J 

where in T g the first factor (2) is due to gluon spin which is twice the quark 
spin, and the second is Casimir scaling factor. 
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One can see that for T c <T < 2T C this ratio is larger than unity, implying 
strong interaction due to magnetic vacuum, especially among gluons. There- 
fore gluons (and quarks) cannot be considered as a gas of weakly interacting 
objects, e.g. for T « T c « 0.2 GeV one has T q «5, ~ 20 and this is 
closer to a liquid rather than a gas. 
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